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Abstract 

We study the decays of top squarks (£1,2) and bottom squarks (61,2) in the Min- 
imal Supersymmetric Standard Model (MSSM) with complex parameters At,Ab, fx 
and M%. We show that including the corresponding phases strongly affects the 
branching ratios of t\ 2 and 612 decays in a large domain of the MSSM parameter 
space. This could have an important impact on the search for t\ 2 and and the 
determination of the underlying MSSM parameters at future colliders. 
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1 Introduction 



Many phenomenological studies on SUSY particle searches have been performed in the 
Minimal Supersymmetric Standard Model (MSSM) with real SUSY parameters. In gen- 
eral, however, some of the SUSY parameters may be complex, in particular the higgsino 
mass parameter /i, the gaugino mass parameters M li2j3 and the trilinear scalar coupling 
parameters Af of the sfermions /. The SU(2) gaugino mass parameter M 2 can be chosen 
real after an appropriate redefinition of the fields. Not only the CP-violating observables 
(such as fermion EDMs) but also the CP-conserving observables (such as cross sections 
and decay branching ratios) depend on the phases of the complex parameters, because in 
general the mass-eigenvalues and the couplings of the SUSY particles involved are func- 
tions of the underlying complex parameters. For example, the decay branching ratios of 
the staus f 12 and r-sneutrino v T can be quite sensitive to the complex phases of the stau 
and gaugino-higgsino sectors pQ. Therefore, in a complete phenomenological analysis of 
production and decays of the SUSY particles one has to take into account that Af, /i and 
M 13 can be complex. In this article based on [2J E| we study the effects of the phases of 
A t , A b , /i and M 1 on the decay branching ratios of the stops ti )2 and sbottoms &i i2 with 
Qi (#2) being the lighter (heavier) squark. We take into account the explicit CP violation 
in the Higgs sector. 



2 SUSY CP Phase Dependences of Masses, Mixings 
and Couplings 

In the MSSM the squark sector is specified by the mass matrix in the basis (g^, g#) with 
q = t or b 

Ml = ( < a >« \ (1) 



q \ a q m g 




with 



m\ L = M 2 Q + m 2 z cos 2/? (if - e q sin 2 W ) + m 2 , (2) 



»§ B = M 2 {(j b} + m| cos 2(3 e q sin 2 6 W + mj, (3) 

/ (A t -fi* cot 0) m t (q = t) 
q 9 { (A 6 -/i*tan/3)m 6 (q = b) { > 

= \a q m q \ e 1 ^ (-7T < (f q < it). (5) 

Here If is the third component of the weak isospin and e q the electric charge of the quark 
q. Mqjj i) and A tjb are soft SUSY-breaking parameters and tan/5 = v 2 /vi with v± (t> 2 ) 
being the vacuum expectation value of the Higgs field (if 2 ). We take A q (q = t,b) 
and fi as complex parameters: A q = \A q \e lipA i and /i = e 1 ^ with — tt < ^fA q ,^ < 
Diagonalizing the matrix one gets the mass eigenstates q\ and g 2 



mi = Af; 





e l</3 « cos 6 q sin 9 q 
— sin 9 q e~ t(pi i cos 9 q 




(6) 
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with the masses and m q2 (m^ < mq 2 ), and the mixing angle 9q 

m i, 2 = \( m l L + m \ R =f \j ( m i - m L.) 2 + 4 IS m J 2 ); ( 7 ) 

5 = tan-^la^l/^-mf^)) (-tt/2 < 5 - < 0). (8) 

The — qji mixing is large if |m? L — rrv- \ < |a ? m g |, which may be the case in the t sector 
due to the large m t and in the b sector for large tan/3 and \fj,\. 

We assume that the gluino mass m g = M 3 is real. We write the U(l) gaugino mass 
Mi as Mi = \Mi\e t(fl (—it < (pi < n). Inspired by the gaugino mass unification we take 
|Mi| = (5/3) tan 2 9 W M 2 and m g = (a s (m g ) / a 2 )M 2 . In the MSSM Higgs sector with 
explicit CP violation the neutral Higgs mass eigenstates H®, and H® (m H o < m H o < 
m H o) are mixtures of CP-even and CP-odd states. For the radiatively corrected masses 
and mixings of the Higgs bosons we use the formulae of Ref . [3] . 

Possible important decay modes of ti t 2 are: 

ti - tg,tx° i ,bxt,biW+,b l H + (9) 
t 2 - tg,tx° i ,bxt,tiZ°,~b 1>2 W + , ^ fc Q ,& 1|2 tf + . (10) 

Those of &i ; 2 are analogous. The CP phase dependences of the decay widths stems from 
those of the involved mass-eigenvalues, mixings and couplings among the interaction- 
eigenfields. In E] it is found that the masses, mixings and couplings of the involved 
SUSY particles (stops, sbottoms, charginos xt ( m xf < m xf^ neutralinos Xj ( m x? < ••■ < 
m^o), and Higgs bosons) can be very sensitive to the CP phases in a large region of the 
MSSM parameter space. For the g« (q = t, b) sectors we find: 

1. The mass-eigenvalues m qi2 are sensitive to the phases (pA q ,f^) via cos(<^4 + f^) 
if and only if |a 9 m 9 | ~ (m| £ + m| H )/2 and \A q \ ~ \^\C q (with C t = cot/3 and 
C b = tan/5) (see Eqs. © , © ) . 

2. The g-mixing angle 9q (given by tan 29 q = 2\a q m q \/(m^ L — m~ R ) ) is sensitive to 
{(PA q ,Vn) via cos(p Aq + ¥V) if and only if 2\a q m q \ > \rri~ L - mfj and \A q \ ~ \n\C q 
(see Eqs. ©,©, ©). 

3. The g-mixing phase tpq in Eq.© is sensitive to (<pA q , w) 5 VA q and y? M if \A q \ ~ \n\C q , 
\A q \ S> \fi\C q and |^4 9 | <C \/J>\C q , respectively. For large squark mixing the term 
oc sin 29 g cos (fq can result in a large phase dependence of the decay widths |2] (see 
Eq.®). 

Therefore, we expect that the widths (and hence the branching ratios) of the ti t2 and bi >2 
decays are sensitive to the phases (<fA tb , </V, (pi) in a large region of the MSSM parameter 
space. 
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3 Numerical Results 



In order to improve the convergence of the perturbative expansion [Sj we calculate the tree- 
level widths by using the corresponding tree-level couplings defined in terms of "effective" 
MSSM running quark masses (i.e. those defined in terms of the effective running 

Yukawa couplings h r t u b n oc m r t ^ 1 ). For the kinematics, e.g., for the phase space factor we 
use the on-shell masses obtained by using the on-shell (pole) quark masses M tjb . We 
take M t = 175 GeV, M b = 5 GeV, m r t un = 150 GeV, and m r b un = 3 GeV. We fix 
\A t \ = \A b \ = \A\ and M 2 = 300 GeV, i.e. m- g = 820 GeV. In our numerical study 
we take tan/3, m^, mj 2 , m^, \A\, (fA t , <fA b , (fn, (pi and m H + as input parameters, 
where 2 and mi are the on-shell squark masses. Note that for a given set of the 
input parameters we obtain two solutions for (Mq, Mfj) corresponding to the two cases 
m t L > m t R and vai h < m; from Eqs. (JU)-(jU) and (J7J) with m t replaced by M t . In the 
plots we impose the conditions as described in El in order to respect experimental and 
theoretical constraints, such as the LEP limits and h — > 57 data. 

The experimental limits on the EDMs of electron, neutron and 199 Hg strongly constrain 
the SUSY CP phases. Here we adopt the scenario of [6 3 where the first two generations 
of the sfermions are very heavy and hence <fA t b , and tp\ are practically unconstrained. 
We have also checked that the electron and neutron EDM constraints at two-loop level 
[7j are fulfilled in the numerical examples studied in this article. 

In Fig.l we show the ipA t dependence of the ti decay branching ratios for (f^ = and 
tt/2 with tan/3=8, (m fl ,m t ~ 2 ,m^) = (400,700,200) GeV, \A\ = 800 GeV, = 500 GeV, 
ipA b =( Pi = ®, and m H + = 600 GeV in the case m t ~ L > m^ R . The case < mi R leads to 
similar results. We see that the t 2 decay branching ratios are very sensitive to (fA t and 
depend significantly on (p^. For large tan/5( > 15) we have obtained results similar to 
those for tan/3=8 [21 El- For ti and b X 2 decays we have obtained results similar to those 
for the t 2 decays (21 El- We have also found that for small tan/3( < 8) the t\ decay can 
be fairly sensitive to tp\ [2]. 

4 Parameter Determination 

We now study to which extent one can extract the underlying MSSM parameters (such as 
tan/3, Atfi etc.) from measured observables (such as masses, cross sections and branching 
ratios). The observables are functions of the MSSM parameters. Therefore, in case the 
number of measured observables is larger than that of the MSSM parameters, the whole 
data set of the observables over-constrains the MSSM parameters. In this case the MSSM 
parameters and their errors can be determined by a x 2 fit to the experimental data of the 
observables. In the fit the errors of the observables are propagated to those of the MSSM 
parameters. The assumptions for the expected experimental errors of the observables 
measured at TESLA, CLIC and LHC are described in Our strategy for the parameter 
determination is as follows: 
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Figure 1: (p At dependence of the decay branching ratios for (p^ = (a) and 7r/2 (b) with 
tan 0=8, (m t - 1 ,m t - 2 ,m Si )=(400,700,200) GeV, \A\ = 800 GeV, |/x| = 500 GeV, VAb = Vl =Q, 
and rriH+ = 600 GeV in the case mjr > rrif . Note that the tiH^s and b\H + modes are 
kinematically forbidden here. 



1. Take a specific set of values of the underlying MSSM parameters. 

2. Calculate the masses of fj, fy, xp Xki H®, the production cross sections for e + e~ — > 

tjtj, and e + e~ — > and the branching ratios of the U and hi decays, and estimate 
the expected experimental errors of these observables. 

3. Regard these calculated values as real experimental data with definite errors. 

4. Determine the underlying MSSM parameters and their errors from the "experimen- 
tal data" by a x 2 fit- 



We consider two scenarios, one with small tan (3 and one with large tan (3. The small 
tan/? scenario is characterized by the following underlying MSSM parameters: Mg, = 
169.6 GeV, M = 408.8 GeV, Mq = 623.0 GeV, \A t \ = \A b \ = 800 GeV, (p At = <p Ab = tt/4, 
Vl = 0, M 2 = 300 GeV, y. = -350 GeV, tan (3 = 6, m- g = 1000 GeV, and m H ± = 900 GeV. 
(Here we do not assume the unification relation between m- g and M2.) The large tan/3 
scenario is specified by: M D = 360.0 GeV, M = 198.2 GeV, Mq = 691.9 GeV, \A t \ = 600 
GeV, ip At = vr/4, \A b \ = 1000 GeV, ip Ab = 3vr/2, ip 1 = 0, M 2 = 200 GeV, /i = -350 GeV, 
tan/5 = 30, m g = 1000 GeV, and = 350 GeV. The resulting values of the observables 
and their expected experimental errors are shown for the two scenarios in We regard 
these calculated values as real experimental data with definite errors. We determine the 
underlying MSSM parameters and their errors from the "experimental data" on these 
observables by a x 2 fit- The results obtained are shown in Table 2 of As one can see, 
all parameters except A b can be determined rather precisely, tan j3 can be determined with 
an error of about 3% in both scenarios. The relative error of the squark mass parameters 
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squared M? ^ ^ is in the range of 1% to 2%. Re(A t ) and Jm(v4 t ) can be measured within 
an error of 2 - 3% independently of tan (3. The situation for is considerably worse: in 
case of small tan (3 one gets only an order of magnitude estimate. The reason is that both 
the bottom squark mixing angle and the bottom squark couplings depend only weakly on 
Ab for small tan (3. In case of large tan f3 the situation improves somewhat in particular 
for the imaginary part of Af,. A x 2 fit using only real MSSM parameters would result 
in a totally wrong parameter determination with a significantly larger value of \ 2 . We 
have found that the analogous fit procedure using only real MSSM parameters gives a 
much larger value for x 2 '- ^X 2 = 286.6 with DOF=61 for the scenario with tan/? = 6 
and Ax 2 = 22.5 with DOF=61 for the scenario with tan/? = 30, where DOF is 'degree of 
freedom'. 



5 Conclusion 

We have studied the decays of ti j2 and &i i2 in the MSSM with complex parameters A tjb , fi 
and Mi. We have shown that including the CP phases strongly affects the branching 
ratios of ti i2 and &i )2 decays in a large domain of the MSSM parameter space. This could 
have an important impact on the search for ti j2 and &i j2 and the determination of the 
underlying MSSM parameters at future colliders. 
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